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Abstract 

The spin-dependent cross sections for semi-inclusive lepton-nucleon scattering are derived in the 
framework of collinear factorization, including the effects of masses of the target and produced 
hadron at finite momentum transfer squared Q^. At leading order the cross sections factorize 
into products of parton distribution and fragmentation functions evaluated in terms of new, mass- 
dependent scaling variables. The size of the hadron mass corrections is estimated at kinematics 
relevant for future semi-inclusive deep-inelastic scattering experiments. 
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I. INTRODUCTION 


Determining the detailed flavor and spin structure of the nucleon remains a central chal¬ 
lenge for hadronic physics into the 21st century. Considerable progress has been made over 
the past two decades in understanding the characteristics of the momentum and spin dis¬ 
tributions of quarks and gluons (or partons) through precise measurements of the nucleon’s 
parton distribution functions (PDFs) in various hard scattering reactions [l, 2|. In addition 
to the traditional inclusive deep-inelastic scattering (DIS), Drell-Yan and other high-energy 
scattering processes, an increasingly important role in this quest has been played by semi- 
inclusive production of hadrons in lepton-nucleon scattering. 

Identification of specific hadrons, such as pions or kaons, in the current fragmentation re¬ 
gion of a deep-inelastic collision serves as a tag of individual quark flavors, which in inclusive 
DIS are summed over. Important insights have been provided through semi-inclusive deep- 
inelastic scattering (SIDIS) experiments on phenomena such as the SU(2) flavor asymmetry 
in the proton sea |3| and the ratio of strange to nonstrange quark distributions |^. From 
experiments with polarized targets, SIDIS data have also provided fascinating glimpses of 


the possible flavor asymmetry in the polarized light-antiquark sea j^, while kaon produc¬ 
tion data has fueled the recent controversy concerning the sign of the polarized strange sea 
g, l7|]. Furthermore, detection of forward baryons (in the center of mass frame) in the target 
fragmentation region of SIDIS is a potentially important avenue for extracting information 
on the pion cloud of the nucleon or the structure of the virtual pion itself 

In more recent developments, detection of both the longitudinal and transverse momen¬ 
tum distributions of hadrons produced in SIDIS measurements of various single- and double¬ 
spin asymmetries has opened u p t he largely unexplored realm of transverse momentum 


dependent parton distributions 12l-ll4|. These reveal an even richer landscape of three¬ 


dimensional momentum and spin distribution of partons in the nucleon, that will be the 
subject of increasingly greater attention at facilities such as Jefferson Lab 1^ and COM¬ 
PASS 16 L and a central component of the science program at the proposed Electron-Ion 


Collider 


17 |. 


The unambiguous interpretation of any SIDIS experiment in terms of leading twist PDFs 
or transverse momentum distributions requires control of various subleading 1/Q^ correc¬ 
tions, such as target mass and higher twist effects, as well as knowledge of the fragmentation 
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functions describing the parton hadronization. For inclusive DIS, the £nite-Q^ corrections 
are known to become important at low values, particularly when the parton momentum 
fraction x is large 18l-l20|. Their effects on global hts of spin-averaged PDFs have been 
systematically studied in recent analyses by the CTEQ-Jefferson Lab (CJ) 


M 




ABM 


25| and JR groups [2^ , and in spin-dependent PDF analyses by 


27 | (and to some extent also by the LSS 28|, BB 29| and NNPDF 


he JAM collaboration 


301 1 groups). 


Typically, the effects of target mass corrections (TMCs) can be computed within a specific 
framework, while higher twist effects, which involve more complicated multi-parton correla¬ 
tions, are parametrized phenomenologically. The standard approach for computing TMCs 
has traditionally been within the operator product expansion, in which the mass corrections 
to inclusive DIS structure functions arise from twist-two quark bilinear operators with an 

n n 

arbitrary number of derivative insertions |18l. l31| . Extending this framework to processes 
involving particles in the final state is problematic, however, which has in practice limited 
the study of hadron mass corrections in SIDIS. 

An alternative framework for TMCs was developed using techniques based on collinear 
factorization (CF) 32h34| . in which the hard scattering is formulated in momentum space 


directly. The method has 
both in unpolarized |35M37 


Deen applied to the computation of TMCs in inclusive scattering, 
and polarized 


DIS, and in semi-inclusive hadron production 
in electron-proton anhilation 39|, 1^. For semi-inclusive hadron production in lepton-proton 
collisions, in contrast to inclusive DIS and e’''e“ anihilation, finite-Q^ corrections can arise 
from both the effects of the target mass and the mass of the produced hadron. While 
earlier analyses 4l|, l42|] considered some of these corrections within the CF framework, the 
phenomenology of the combined effects of the target and produced hadron masses - which 
we refer collectively as “hadron mass corrections” (HMCs) - was systematically explored in 


Ref. 


43j for unpolarized scattering. 


In this work we extend the analysis of HMCs to the case of spin-dependent SIDIS at 
hnite Q^. Because high energy spin-dependent data are generally more scarce than spin- 
averaged cross sections, a significantly larger fraction of the world’s data set used to constrain 
spin-dependent PDFs lies in the low-Q^ region {Q‘^ ~ 1 — 2 GeV^). While target mass 


corrections have been incorporated in some global spin-PDF analyses 


27 


analyses which have included polarized SIDIS data 


28|, 


30l |. none of the 


44l | have accounted for HMCs. With 


the increasing precision of new polarized measurements, and the consequently more accurate 
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determination of spin-dependent PDFs, it is imperative to reliably account for subleading 
corrections which could impact the extraction of the leading twist distributions. 

In Sec. mi we outline the formalism used to compute the SIDIS cross sections at finite val¬ 
ues of within the collinear approximation in the presence of target and produced hadron 
masses. For completeness, we consider both polarized and unpolarized scattering, since the 
latter enters the calculation of the measured polarization asymmetries. In Sec. IIIBI we re¬ 
view the collinear formalism and its application to hadron production in SIDIS. Expanding 
the hadronic tensor in terms of quark correlators, in Sec. Ill Cl we derive semi-inclusive cross 
sections, which at leading order are given by factorized products of PDFs and fragmenta¬ 
tion functions expressed as functions of modihed scaling variables. The relative importance 
of the HMCs is explored numerically in Sec. IIIIl where we quantify the dependence of the 
hnite-Q^ cross sections on the kinematical variables and estimate the corrections for specific 
current and future experiments. Finally, we conclude by summarizing our results in Sec. IIVI 


II. SEMI-INCLUSIVE SCATTERING WITH MASS CORRECTIONS 

The semi-inclusive lepton-nucleon scattering process is illustrated in Fig. [H where the 
incident lepton (with momentum i) scatters from an initial state nucleon (p) to a recoil 
lepton {£') via the exchange of a virtual photon (g), producing a final state hadron h (with 
momentum ph). In this section we first review the external kinematics and choice of variables, 
before outlining the collinear factorization framework for describing the hard scattering 
process. After defining the hadronic tensor in terms of quark-nucleon and quark-hadron 
correlation functions, we derive expressions for the spin-averaged and spin-dependent cross 
sections in terms of parton distribution and fragmentation functions at leading order in the 
strong coupling constant and for hnite values of Q^. 


External kinematics 


We expand the four-momenta of the external particles in terms of coplanar light-cone 
unit vectors n and n, satisfying = 0 and n -n = 1 3^. The “plus” and “minus” 

components of a four-vector are dehned hj v~^ = v ■ n = (v^ + v^)/\/2 and v~ = v -n = 
(u° — v^) / y/2. We work in a class of reference frames where the initial nucleon and virtual 
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FIG. 1; Semi-inclusive deep-inelastic lepton-nucleon scattering with production of a final state 
hadron h. The external momenta of the incident (i) and scattered {£') leptons, virtual photon 
(q), target nucleon (p) and observed hadron (ph) are labeled explicitly. The unobserved final state 
hadrons are labelled by X. 


photon momenta are coplanar, so that 


+ 


2p+ 




+ 




2^p+ 


■ 


( 1 ) 

( 2 ) 


where M is the nucleon mass, and the scaling variable ^ The 

(p, q) collinear frames are related to each other by a boost of parameter p"*" and contain, 
in particular, the target rest frame (where the nucleon p"*" component is p"*" = M/\/2) and 


the Breit 
of Ref. 

45 


31 


rame (p’*' = Q/{y/2^)). For other choices of reference frames see Appendix A 


43|. In the collinear (p, g) frame, ^ is identical to the Nachtmann scaling variable 

m 




2x 


B 


( 3 ) 


1 + v/1 + 4x|MVQ2’ 

where xb = Q^/2p • g is the Bjorken scaling variable, and in the Bjorken limit {Q‘^ and 
g“ —)• oo with Xb hxed) one hnds that ^ Momentum and baryon number conservation 

(and, for K production, strangeness conservation) impose an upper limit on x^, xb < 
where 

mhirrih + 2 Mb) + 


= 1 + 




( 4 ) 


B V 

with Mb = M for tt production and Mb = Ma (the lightest hyperon mass) for K production. 
This value of x^ corresponds to a hnal state consisting of the nucleon and hadron h at rest 
in the hadron’s rest frame. 
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For spin-dependent scattering, the polarization vector of the initial state nucleon can 
be parametrized as 




M\ 


p. 




■ 


Cti 

I ) 


(5) 


and satishes the conditions p ■ S = 0 and = — 1. In the case of a longitudinally polarized 
initial state nucleon (which we consider in this work), one has Sl = ±1 and the transverse 
spin vector = 0 . 

The incident and scattered lepton momenta can be decomposed as 

= VP'^n^ + (l+ (6) 

(7) 


iJ2ip+ 

= {p + ^)p'^n'^ + T ^ I 




^ 2^p+ 

where is the lepton transverse momentum four-vector, p = /p'^ is the lepton momentum 

fraction, and we assume massless leptons, = £'‘^ = 0. After some algebra one can show 
that 


p = 




22/(1 + 7 ^) 


(2 - 2 /)Vl + 72 - y{l + 7 ^) 


( 8 ) 


where y = q ■ p/i ■ p and 7 ^ = In the target rest frame y = u/E is the fractional 

energy transfer from the lepton to the target, with u and E the virtual photon and incident 
lepton energies, respectively. The magnitude of the lepton transverse momentum is set by 
four-momentum conservation. 


p 




= -4 = -7 1 + 7 Q • 


p 


For the hadron produced in the hnal state, the momentum is parametrized as 


(9) 


( 10 ) 


where C,h = Ph/Q~ is the scaling fragmentation variable, and the hadron transverse momen¬ 
tum four-vector satishes phi. ■ n = phi_ • n = 0, with norm p^j_ = —Ph±- The squared 

transverse mass of the produced hadron h is dehned by m^_|_ = -|- Ph±, where rrih is the 

mass of the hadron. The variable (h can be related to the fragmentation invariant 


Ph-P ^ , 

Zh = -= u + 

q-p ^ 


( 11 ) 


a 

In the target rest frame it coincides with the ratio of the produced hadron to virtual photon 
energies, Zh = E^/v^ which is frequently used in experimental analysis of SIDIS data. In the 
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Bjorken limit the fragmentation variable (h Zh, while at hnite one has 

Zh i i ^ , /-, 4:x\M‘^rrii 


Ch = 


2 xb 


1 + W1 — 






( 12 ) 


Since the produced hadron’s energy is bounded from below by Eh > mhx, one can show 
that Zh > -2™™, where 


= 2xb- 


Mmh 




(13) 


Combining Eqs. (IT^ and flT^ . one can show that the corresponding minimum value of C,h is 
given by = ^Mruh/Q"^■ In the target rest frame, 2 ;™’° corresponds to the hadron h pro¬ 
duced at rest, with the remaining hnal state hadrons moving collectively in the direction of 
the virtual photon. At the other extreme, conservation of four-momentum, baryon number, 
and (for K production) strangeness impose the upper limit Zh < where 

M(Mb - M) 


= 1 - 2xb- 




(14) 


with again = M for h = n and = M\ for h = K. This limit corresponds to diffractive 
production of the observed hadron with maximal energy. As —)■ cx), both the upper and 
lower limits become independent of 2 ;™™ —)■ 0 and z™^^ —)■ 1. 

As an alternative to the fragmentation invariant Zh, one can define the invariant momen¬ 
tum fraction 


_‘2ph-q _ ^ 
Zp. — o — Ch 


m 


h± 


r O.Q=’ 

which is used in the study of hadron production in e+e“ collisions. The choice of this variable 
avoids mixing inclusive {xb) and semi-inclusive {zh) variables as in Eq. flTT]) . It also allows 
a clean separation of the current {zp > 0) and target {zp < 0) fragmentation regions in the 
Breit frame, in which = Ph/q^ is ratio of the longitudinal components of the hadron 
and photon momenta. In the current region, where the observed hadrons are produced with 
longitudinal momentum in the direction of the virtual photon, Ch can be written in terms 
of the variable as _ 

L 4:ml 


Ch = 


1 + \ 1 + 


h± 


( 0 ) _ 


(16) 


Hadrons produced in the current region have Ch > Ch = Ch{zp = 0), where 

.(0) _ 

~ Q ' 


( 17 ) 
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FIG. 2: Finite-Q^ fragmentation variable Ch versus Zh for the semi-inclusive production of (a) pions, 
h = TT and (b) kaons, h = K, at fixed values of xb = 0.3 (blue curves) and 0.6 (red curves) for 

= 1 (solid curves) and 5 GeV^ (dashed curves). The curves are shown only in the kinematically 
allowed Zh regions, and the boundaries between the current {(h > Ci^^) target {Ch < Ch'^) 
fragmentation regions are indicated by the open circles. 

Note that which reflects the fact that a hadron produced at rest in the target 

rest frame belongs to the target region. Finally, in the Bjorken limit all three fragmentation 
variables become equivalent, (k —t Zh —t Ze, and the current region extends down to the 
smallest values of z^, —>■ 0. 

The relation between the variable and z^ is illustrated in Fig.[2]for several hxed values of 
xb and Q^. At hnite the kinematically allowed regions of Zh are determined by Eqs. flT3|) 
and (na, and the boundaries between the current and target fragmentation regions occur 
at Ch = For the production of pions, at low xb ^ 0.3 and high > 5 GeV^ the 
differences between the two variables are almost negligible, and begin to be noticeable only 
for the highest Zh values at lower or higher xb- At these kinematics, pions are produced 
in the current region for Zh ^ 0.08 — 0.15. At high xb ~ 0.6 and low Q'^ ^ 1 GeV^, the 
effects are more pronounced, with deviations of ~ 30% as —)■ 1, and current fragmentation 
begins at a slightly higher Zh- The effects of the kinematic lower limit in Zh [Eq. (1T3|) ] is 
noticeable only at low and high xb- 

For kaons, the effects at high and low xb are again negligible, although the larger K 
mass enhances the differences relative to the pion at the same kinematics. In particular, at 
low = 1 GeV^ the lower limit on Ch is dramatically increased, and the current region is 
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FIG. 3: Semi-inclusive deep-inelastic lepton-nucleon scattering with production of a final state 
hadron h at leading order in as- The internal momenta of the initial (k) and scattered quarks {k') 
are labeled explcitly. The intermediate state X' represents a nucleon with a quark removed, and 
the state Y' results from the fragmenting parton with the hadron h removed. The dashed vertical 
line represents the cut in the forward scattering amplitude. 

pushed to higher values of Zh- At sufficiently large xb, the phase space for K production 
eventually vanishes; for = 1 GeV^, for example, no kaons can be produced in the current 
fragmentation region with xb ^ 0.66, and no kaons can be produced at all for xb ^ 0.8. 


B. Collinear factorization 

At leading order in the strong coupling Os, the SIDIS reaction proceeds through the 
hard scattering of the virtual photon from an initial state quark with momentum /c to a 
quark with momentum k' = k + q, which then fragments to a hadron h, as illustrated in 
Fig. |3l Higher order processes involving gluon radiation and scattering from qq pairs can be 
considered, but for clarity of the derivation of the finite-Q^ corrections we restrict ourselves 
to the leading order calculation. 

The parton four-momenta can be parametrized, in analogy with the external variables in 
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Sec-EAl in terms of the light-cone vectors n and n as 


^2 

= xp'^ - ^n^ + kk, (18) 

2xp+ 

yt. ^ ^ ^ ^ glL (19) 

2ph /z z 

where x = k~^/p'^ is the light-cone momentum fraction of the nucleon carried by the struck 
quark, and z = p'^/k'~ is the light-cone fraction of the fragmenting quark carried by the 
hadron h. The parton transverse momentum four-vectors k± and k'^ are orthogonal to n 
and n, k± ■ n = k± ■ n = {) and k'^ - n = k'^ - n = with norms k\ = —k\ and k']_ = 

The essence of the collinear factorization approach is to expand the hard scattering term 
around “collinear” momenta k and k', 


k^ 


k'^ 


Xp+ _j_ 


2xp+ 




+ Pl±h^ 


z z ■ 


( 20 ) 

( 21 ) 


where the initial and hnal parton virtualities k^ and k'"^ are kept for generality. In this 
approximation the transverse momentum of the initial parton is neglected, and the transverse 
momentum of the hnal parton is taken along the direction of the hadron h, making these 
collinear to the proton target and produced hadron, respectively. 

Dehning the partonic analog of the Bjorken scaling variable at hnite by 

e 1 


X = 


Q 


2k ■ q XI — ^ 2 /^ 2 ^ 2 ; 2 Q 2 
one can show that x lies within the range Xmin < x < fmax, where [37|, 143 
1 1 2Mmh + k'^ 1 


Xr, 


1 

Xb 




Xr, 


= 1 + 


ml 


^(1 




xC,hQ‘ 


( 22 ) 


(23) 


Here the lower limit Xmin corresponds to the minimal mass of collinear spectator partons, 
while the upper limit Xmax arises from the minimum value of the current jet invariant mass. 
These limits are consistent with the limit on xb in Eq. (j4]) for any > x{C,h — ^)Q'^/^i nnd in 
the Bjorken limit the range of x reduces to being between xb and 1, as required. Analogous 
limits can be found for the partonic fragmentation variable z, Zmin ^ z < 1, where 

^min C/i ^ X 
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For the practical implementation of the hnite-Q^ kinematical limits, a choice of the ini¬ 
tial and scattered parton virtualities /c^ and /c'^ must be made. For the initial parton the 
requirement that the collinear parton virtuality is independent of x leads to the restriction 
/c^ > 0. For light, bound initial-state partons (/c^ < 0) this constrains therefore the collinear 
expansion to be around = 0. Determining the virtuality of the scattered parton is 
generally less clear, on the other hand. For the leading order hard scattering in Fig. [3l 
conservation of four-momentum and the condition = 0 constrain the parton momentum 
fraction x to equal where 


a = e 1 + 




(25) 


m 


In order for x to respect the limits in Eq. fl2^ . the choice k'"^ = rn‘l/(h was proposed 
Ref. 1 ^, giving = ^(1 + m\/(hQ^). Larger k'^ values would also allow x to fall within 
the bounds in Eq. fl23|) : however, the above choice is the closest to the physical quark mass. 


and is the one adopted in our numerical analysis here. The dependence of the calculated 
cross sections on the choice of k'^ is explored further in Sec. IIIII 


C. SIDIS at finite Q‘^ 


The cross section for the SIDIS process can be written as a product of hadronic (W^'^ 
and leptonic {L^^) tensors [l3|, 

da 




dxBdyd^Ph/2Eh 


(26) 


where a is the electromagnetic hne structure constant. The leptonic tensor can be computed 
straightforwardly from QED, 


f, A) = 2i/^ + 2iJ'^ + + 2t\ (27) 

where A is the lepton helicity. 

The semi-inclusive hadronic tensor is defined in terms of matrix elements of the electro¬ 
magnetic current operator between the initial state nucleon with spin S and the hnal 
state with a hadron h and unobserved hadrons X, 

2MW>^''{p,S,q,ph) = J + Q - Px - Ph) 

X {N, 5| J'^(0)|h, X) {h,X\r{0)\N, S), (28) 
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where we use the shorthand notation cPpx/‘2Ex = Yliex with px the total mo¬ 

mentum of teh unobserved hadrons. The hadron tensor can be expressed at leading order 
in as in terms of quark-quark correlators <f>jand A^, associated with the quark distribution 
and fragmentation functions, respectively 


33 




m 


2MW^'^{p,S,q,pu) = Y,el / d^k d^k' 5^^\k + q-k')i:T S,k) ^ ,ph) Y] , (29) 


where the sum is taken over quark and antiquark flavors q. Note that since the parton 
momenta in the (5-function have been approximated by their collinear components, the in¬ 
tegrations over dk~ d?k\_ and dk'^ act directly on the correlators and . 

The correlator associated with the parton distribution function is defined, in the light- 
cone gauge, as 




(2ir) 


X' 


2Ex' 


where il^q is the quark held operator, and Ex' and px' are the energy and momenta of the 
intermediate state X' corresponding to a nucleon with a quark removed. Similarly, for the 
quark fragmentation correlator one has 




(27r) 


Y' 


where Ey and py are the energy and momenta of the state Y' resulting from the fragmenting 
quark with the hadron h removed. 

The leading twist contributions to the SIDIS cross section can be extracted by expanding 
the integrated correlators and selecting the terms linear in the light-cone vectors. For the 
parton distribution function, one hnds 

<hq(x) = J dk~d‘^k±^q{p,S,k) = ^q{x)^ + Ag(x) 7 ® ^..., (32) 

where the hrst and second terms dehne the spin-averaged, q{x), and spin-dependent, Aq{x), 


distribution functions, and the ellipsis indicates higher twist contributions 
mentation function D^[z) is analogously dehned from the A^ correlator by 

x,{z) ^ i / dk’+d'XX,(k',p,) = + ... 


131. The frag- 


(33) 
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Inverting Eqs. (I32|) - (l33|) . one can write the qnark distribntion and fragmentation functions 
explicitly by projecting with the appropriate Dirac operators, 


q{x) = 


Aq{x) = 


2 J 
1 

47r 

D 

1 

Att 

Stt 


=xp^ 


I dk d^fcxTr [7+$g(p, S', 

j dw- 7+ V>,(tt’"n)|iV>, 

( dk-d?k^ Tr [7'7*4,(p, S, k)] 

Jdw~ (Af|^g(0)7^7"'"V'g(w^p)|Af), 

f dk'*d^k'^Tv [7-Aj(fc'.p^)]^,_____^^ 

j dw* e‘^<‘kl‘>\0\i,,(w*n)\hX){h,Y'\i’,(li)r\0), 


(34) 


(35) 


(36) 


Y' 


where are light-cone coordinates. (For ease of notation we omit the dependence in 
the arguments of the quark distribution and fragmentation functions.) The fragmentation 
function D^{z) here is dehned with the standard normalization, lo ^ ~ ^ — I- 

To compute the hadronic tensor in Eq. ([29]), we can decompose the + q — k') func¬ 
tion along the +, — and transverse components of the momenta. The ^-functions for the -|- 
and — components constrain the partonic variables to x = and z = (h, respectively. The 
(5-function for the transverse component forces the transverse momentum of the produced 
hadron h to vanish, ph± = zk'j_ = 0. Nonzero transverse momentum hadrons can be pro¬ 
duced via higher order perturbative QCD processes, or from intrinsic transverse momentum 
in the parton distribution functions themselves 13| . With these constraints, the hadron ten¬ 
sor can then be factorized into products of parton distribution and fragmentation functions 
evaluated at and Ch, respectively, 

Ch 


2MW^‘''(p,S,q,p,) = y^ep®(px)('IV[^7W9({/.) 


+ SilV A«(&) + ■ ■ •)£>;(().)■ (37) 


The main effect of the hadron masses at hnite kinematics is therefore a replacement of 
the Bjorken limit scaling variables xb and Zh by their hnite-Q^ analogs. However, since 
depends explicitly on m/i, and depends on Zh and xb, the scattering and fragmentation 
parts of the hadronic tensor at hnite are not independent. 

Contracting the hadronic tensor in Eq. fl37D with the leptonic tensor in Eq. fl27|) enables 
the leading order spin-averaged (ah) and spin-dependent (Aah) SIDIS cross sections to be 


13 




written in terms of parton distribntions evalnated at the new scaling variables, 


1 


= t: 




AcTfc = 


2dxB dQ"^ dzh 
dal^-^^ 


~^h'i 


27r«^ y'^ 

i- e 
4:71 


^(yh, 


dxsdQ'^dzh \ — e 

where the rednced nnpolarized and polarized cross sections are dehned as 


(38a) 

(38b) 


A<ti = A5^e5A,(&,Q2)Dj(a,Q"). 


(39a) 

(39b) 


In Eqs. fl38l) the arrows denote the spins of the lepton and target nucleon, and the depen¬ 
dence of the functions on the scale is made explicit. In Eqs. (|38|) the variable 

1 - y - vW!^ 


e = 


1 -!/ + !/^(l+ p^)/2 


(40) 


is the ratio of longitudinal to transverse photon flux, and in Eqs. fl5U]) is a scale dependent 
Jacobian factor, Jh = dC,h/dzh = (1 — Af^^^/(5^)/(l — , with ^ 1 at large 

Q2. 

Note that at the maximum value of allowed for SIDIS [see Eq. (jl])] the hnite-Q^ 
variable satishes = x^^^) < 1. As in the case of inclusive DIS |37|, the 

SIDIS cross section therefore does not vanish as xb —)■ which reflects the well-known 


threshold problem in which the leading twist structure function is nonzero for > 1 


48- 


51| . Analogously, for the £nite-Q^ fragmentation varialble one has C/i < Ch{,Zh = < 1, 

and since the fragmentation function does not vanish as Zh —t the perturbatively 

calculated SIDIS cross section can also exceed the fragmentation threshold. 

Before exploring the dependence of the SIDIS cross sections on the £nite-Q^ scaling 
variables in the next section, it is useful to hrst establish the purely kinematic corrections, 
independent of the PDFs and fragmentation functions, which augment the hnite-Q^ results 
from their scaling limit. In Fig. |4]we show the Zh dependence of the Jacobian factor in 
Eqs. fl3^ at hxed values of xb and for the case of pion production {h = tt). For xb = 0.3 
the factor deviates very little for unity over most of the range of Zh-, with an upturn only at 
small Zh, Zh < 0.1 for > 1 GeV^. At higher xb values the hnite-Q^ effects are more visible, 
with Jh spiking above unity at Zh < 0.2 for xb = 0.8 and = 1 GeV^, and decreasing to 
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FIG. 4: Jacobian factor Jh versus Zh for pion production {h = vr) at fixed values oi xb = 0.3 (blue 
curves) and 0.8 (red curves), at = 1 (solid curves) and 5 GeV^ (dashed curves). The open 
circles denote the boundary between the target (small z^) and current (large Zh) fragmentation 
regions. 


~ 25% below unity at large Zh- This behavior should be kept in mind when assessing the 
numerical effects of the HMCs in the cross sections in the next section. One should also note 
that the upturn in is almost entirely conhned to the target fragmentation region (small 
Zh), where the validity of calculations based on the perturbative handbag diagram in Fig. [3] 


is more questionable, and factorization in terms of fracture functions 
appropriate. 


52 


54| may be more 


III. PHENOMENOLOGICAL IMPLICATIONS 

In this section we examine numerically the phenomenological consequences of the hnite- 
rescaling of the SIDIS cross section derived in Sec. HIl and explore their impact on future 
hadron production experiments. 

A. Kinematical dependence of HMCs 

To disentangle the separate HMC effects in the SIDIS cross sections arising from the 
PDFs and fragmentation functions, in Fig. [5] we illustrate the ratios of PDFs at hnite 
to the corresponding massless limit distributions. For a systematic comparison we consider 
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FIG. 5: Ratios of (a) spin-averaged q = u + d and (b) spin-dependent Aq = Au + Ad isoscalar 
PDFs to the corresponding massless limit distributions, q^^'^ and Aq^^\ as a function of xb at 
various fixed Q^. The finite-Q^ scaling variable here is evaluated for rrih = and = 0.2. 

both the spin-averaged isoscalar q = u + d and spin-dependent Aq = Au + Ad distributions 
at several hxed values of from = 1 GeV^ to 20 GeV^. For the spin-averaged and 


spin-dependent PDFs we use the leading order GT [5^ and LSS [5^ parametrizations, 


respectively, and evaluate the scaling variable for the case of pion production, rrih = 
and for a typical value of (h = 0.2. The results for other (h values are similar, essentially 
given by the PDF evaluated at a rescaled value of [see Eq. fl2^ ]. 

The most dramatic feature in the ratios is the steep rise at large xb, which sets in at 
smaller xb values for decreasing Q^. The results are qualitatively similar for the unpolarized 
and polarized distributions, with the rise delayed to slightly larger xb for the latter. The 
differences between the unpolarized and polarized PDF ratios for the most part reflect the 
differences in the shapes of the respective input PDFs. However, the qualitative features of 
the results do not change when using other leading order distributions, such as the unpolar¬ 
ized GJR and polarized BB 58] parametrizations. Generally, the behavior of the ratios 
observed in Fig. [5] is reminiscent of that found in previous studies of TMGs for inclusive DIS 

Q. 

For the fragmentation functions, ratios of the hnite-Q^ isoscalar functions for 7r+ -|- vr” 
production to those in the massless limit are displayed in Fig. [6] for hxed values of xb 
and Q^, using the leading order HKNS parametrization ^ for the fragmentation functions. 


At Xb = 0.3 the fragmentation function ratio at = 1 GeV^ is enhanced by ~ 20 — 30% 
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FIG. 6: Ratio of vr"*" + 7r“ isoscalar fragmentation functions D to the corresponding massless limit 
functions at fixed values for {a) xb = 0.3 and (h) xb = 0.8. The open circles denote the 
boundary between the target and current fragmentation regions. 


for Zh < 0.7, before rising rapidly as z/i —?■ 1. The effect is less pronounced with increasing 
Q^, with a smaller enhancement of the ratio and a delayed (though even more dramatic) rise 
at large Zh- The fragmentation function ratios at xb = 0.8 in Fig. [6](b), on the other hand, 
display a significantly stronger enhancement, particularly at the lowest value. Here the 
effect is about an order of magnitude larger, and features a striking upturn at Zh < 0.2, 


where the finite-Q^ fragmentation function becomes several times larger than the high-Q^ 


limit. As outlined in Ref. 


43| . this arises from the shape of the fragmentation function at 


finite-Q^ kinematics. The general features of the results, however, remain unchanged if one 
uses the KKP parametrization for instance. 

In particular, expanding the hadron mass corrected fragmentation function in a Taylor 
series about the massless limit, one can write the corrected to uncorrected ratio as 


D{zh) dzh D{zh) 


(41) 


The Zh dependence of the correction is thus determined by the negative shift in the fragmen¬ 
tation variable {Ch — zh) and by the Zh slope of D{zh). Since the pion fragmentation function 
is generally a decreasing function of Zh at small Zh, the ratio is driven upward as Zh 2 :™“, 
where \Ch — Zh\ is maximum. Note that for kaons and protons, in contrast, the slope at 
small Zh can be positive, which would lead to a suppression of the mass corrected function 
at Zh ~ other extreme, in the exclusive production limit the fragmentation 


17 
















FIG. 7: Ratios of spin-averaged {ah/a^\ thin lines) and spin-dependent thick lines) 

cross sections with and without HMCs for semi-inclusive 7r+ -|- tt~ production as a function of Zh, 
at fixed for (a) xb = 0.3 and (h) xb = 0.8. Note the logarithmic scale on the ordinate in (b). 
The open circles denote the boundaries between the target and current fragmentation regions. 

function ratio becomes divergent for the same reason as the PDFs; namely, the functions in 
the scaling limit vanish as xs —?■ 1 or z/j ^ 1, whereas the finite-Q^ scaling variable and the 
corresponding rescaled functions remain hnite. 

Combining the effects of the HMCs in the parton distribution and fragmentation func¬ 
tions, in Fig. Owe show ratios of the SIDIS spin-averaged and spin-dependent cross sections 
with and without HMCs as a function of Zh, for several fixed values of and xb- Specifi¬ 
cally, we consider scattering from a proton target, with the production of 7r+ -|- tt” mesons 
in the final state. The massless limit cross sections and Act® are defined by taking 
the high-Q^ limits of the scaling variables in the arguments of the PDFs and fragmentation 
functions, = ati^h ^ xb, Ch -t Zh) and Acr® = Aah{^h -t xb, Ch -t Zh). 

Overall, the Zh dependence of the cross section ratios follows the trends indicated in 
Figs. [5] and E] for the PDF and fragmentation function ratios, with strong enhancement of 
the finite-Q^ cross sections at large Zh, and decreasing effects at higher Q‘^. At the lower xb 
value {xb = 0.3 in Fig. [71)a)), the HMC effects in the spin-dependent (thick lines) and spin- 
averaged (thin lines) ratios are very similar, which reflects the qualitatively similar shapes of 
the u and Aw distributions at intermediate xb- (Note that the fragmentation functions for 
TT production are the same for the spin-dependent and spin-averaged cross sections.) While 
small differences are visible at Q‘^ = 1 GeV^, at the higher values the unpolarized and 
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FIG. 8: Ratio of the (a) spin-averaged and (b) spin-dependent cross sections with and without 
HMCs for semi-inclusive production of various hadrons h {h = vr'*' -|- vr” or K'^ + K~), for xb = 
0.3 and = b GeV^. The open circles denote the boundary between the target and current 
fragmentation regions. 

polarized ratios are almost indistinguishable. The differences are more striking at larger xb 
{xb = 0.8 in Fig. O^b)), where the effects on the spin-averaged cross section are somewhat 
larger than on the spin-dependent cross section. This stems directly from the delayed rise 
above unity of the PDF ratio in Fig. [5](b) at high xb values compared with the 

corresponding ratio in Fig. [5t^a). 

While most of the existing SIDIS data have involved the production of charged pions, 
the detection of heavier mesons and baryons can provide complementary information on 
the flavor and spin structure of PDFs, as well as on the dynamics of hadronization. The 
production of kaons, for instance, tags strange or antistrange quarks, and has been used 
with polarized targets as an independent means of determining the As distribution in the 
nucleon Q|, and for unpolarized scattering to determine the magnitude of the s distribution 
at small xb jJ]- In Fig. [H] the spin-averaged and spin-dependent cross section ratios with 
and without HMCs are shown for the production of charged pions (7r+ -|- tt”) and kaons 
{K~^ + K~) at Xb = 0.3 and = 5 GeV^. The effects are enhanced signihcantly with 
increasing hadron mass, particularly at low Zh values, mostly because of the (1 -|- m\/(hQ‘^) 
factor in the variable in Eq. ([25]). Increasing values of m\/C,h shift the argument of the 
PDF to higher xb-, where the smaller magnitude of the distributions effectively suppresses 
the cross section. 
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FIG. 9: Ratio of spin-averaged cross sections with and without HMCs for the production of (a) 
pions and (b) kaons, for different choices of the scattered parton invariant mass k'‘^ at = 1 GeV^ 
(thick lines) and = 5 GeV^ (thin lines) for = 0.3. The open circles denote the boundary 
between the target and current fragmentation regions. 


This phenomenon inherently arises from the choice of invariant mass squared for the 
scattered quark, discussed in Sec. Ill B1 While the choice of the mass becomes irrelevant at 
high the dependence on k'^ can be appreciable at low values, as Fig. [9] illustrates. 


Here the spin-averaged cross section ratios computer 
those for massless partons, /c'^ = 0, as used in Ref. 


with /c'^ = m\/C,h are compared with 


40| . For the production of pions, the 


dependence on the quark virtuality is negligible ai = 5 GeV^, but becomes evident at 
lower Q‘^ for small Zh, Zh ^ 0.5. Overall, the cr/i/a® ratio is closer to unity for the preferred 
choice of = fn\/C,h (see Sec. Ill Bll . with greater deviations for the /c'^ = 0 choice. 

For kaon production the effects are expectedly larger, with significant dependence on the 
quark virtuality at = 1 GeV^ for most of the Zh range. Interestingly, the sign of the 
correction is different for the two choices of quark mass for Zh ^ 0.8. By = 5 GeV^ 
the dependence on is weaker, except at Zh ^ 0.3 where visible differences persist. These 
results suggest that care must be taken when extracting PDF information from low-Q^ SIDIS 
data at the extremeties of the Zh spectra, particularly for heavier produced hadrons such 
as kaons. Gaution must also be exercised when including data in the target fragmentation 
region at small Zh, where factorization based on the use of fragmentation functions becomes 
more questionable. In this region utilization of data may require the fracture functions 
formalism as discussed in Refs. 
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FIG. 10: Relative effect of HMCs on the spin-averaged SIDIS cross section as a function of Zh 
or TT p roduction from protons at kinematics typical of future 11 GeV Jefferson Lab experiments 
61I463|. The relative HMG correction is evaluated for each of the three xb values 

{xb = 0.275, 0.575 and 0.725) at two fixed values of (indicated in the legend) from 1.5 GeV^ to 
9 GeV^. 


B. Mass corrections for specific experiments 


The relevance of the HMCs to future SIDIS experiments is illustrated in Fig. [TOl where 
the difference between the hnite-Q^ and Bjorken limit cross sections, 


= a. 


a 


(0) 
h y 


(42) 


is evaluated relative to the 
Jefferson Lab experiments 


inite-Q^ cross section at kinematics typical of planned 11 GeV 


6ll-l63|. The cross sections are computed using the same spin- 


averaged PDF (CT ^) and fragmentation function (HKNS 59|) parametrizations discussed 
in Sec. IIII Al above. 

The effects are pronounced mostly at large Zh, where the ratio of uncorrected to corrected 
cross sections /(Th 0. This is directly correlated with the behavior of the fragmentation 
function ratio in Fig. [6l which diverges as —?■ 1 because at finite the fragmen¬ 

tation variable (h < ^ Sii Zh = 1. The effects at lower Zh are stronger with increasing xb and 
with decreasing Q^. In the range 0.3 ^ Zh ^ 0.6, which is typical for the coverage expected 
in the future experiments, the HMCs are < 10% at xb = 0.275 (where the corresponding 

is between 1.5 and 3.5 GeV^), but increase to ~ 40% at xb = 0.725 (for Q‘^ between 7.5 
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FIG. 11: Effect of HMCs on the (a) unpolarized and (b) polarized SIDIS cross sections for vr 


production from protons at typical 11 GeV Jefferson Lab kinematics 


6Jl |. The cross sections with 


(solid lines) and without (dotted lines) are evaluated at selected xb (and corresponding Q^) values, 
XB = 0.275 (Q2 = 1.5 GeV^), xb = 0.575 (Q^ = 3.5 GeV^), and xb = 0.725 {Q^ = 7.5 GeV^). 

and 9 GeV^). 

The effect of HMCs on the actual (reduced) spin-averaged ah as well as spin-dependent 
Act/i cross sections is shown in Fig. [TTl where the cross sections are calculated at the same 
kinematics as in Fig. [101 The impact of HMCs on the cross sections are more pronounced 
for increasing xb and decreasing Q^. 

For spin-dependent scattering, the effects on the semi-inclusive polarization asymme¬ 
try can also be quantified by defining the difference with respect to the massless limit 
asymmetry 

^(HMC)^h ^ A’l-Af\ ( 43 ) 

where (neglecting the transverse g 2 contribution) A^ = Hy/ZJ. Here = Aah/2ah is the 
parallel asymmetry, and ZJ = (1 — (1 — y)e)/{l + sR) is the photon depolarization factor, 
with R the ratio of longitudinal to transverse photoproduction cross sections. At leading 
order, the polarization asymmetry is then given in terms of the ratios of sums of polarized 
and unpolarized PDFs, 

VT^ E, Ag(e,, Q 2 ) ZJ^"(a, 


A1 = 


(44) 


D E,e2g(a,Q2)ZJ,^(a,Q2) ’ 

where the kinematic prefactor \/l — /D = + 7 ^/(1 -|- 'y‘^y/2). In the Bjorken limit 
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FIG. 12: Effect of HMCs on the SIDIS polarization asymmetry for vr'*' production from protons, 
expressed as (a) relative shift and (b) effect on the asymmetry directly. The asymmetries 


631 ] , for several values of xs and , as 


are evaluated at typical 11 GeV Jefferson Lab kinematics 
in Fig. [TTJ 

this becomes unity, but at finite it represents the projection of the longitudinal lepton 
polarization along the virtual photon direction, cos 9sj^ = D /\/l — where 9sj^ is the angle 
between the lepton and photon spin vectors. 

The corrections to the asymmetry gg actual asymmetry are 

shown in Fig. [12] as a function of Zh for the production of mesons from a proton target, 
at hxed xb and values corresponding to planned 11 GeV Jefferson Lab experiments 631. 
At low Xb values the differences are very small except at very small Zh^ where the effects 
increase as z/i ^ 0. As for the spin-averaged cross sections in Fig. [TOl the effects increase 
with increasing xb and with decreasing Q"^. At the highest xb value, xb = 0.725, the 
asymmetry is reduced by ~ 0.03 for Zh ^ 0.3. If the asymmetry were to approach unity 
in the xb —)■ 1 limit, this would amount to a correction of ~ 3% — 4%. Somewhat larger 
corrections are obtained using the GJR and BB 58] PDF parametrizations, although for 
this combination the asymmetry is not guaranteed to respect the positivity constraint 
at large xb- On the other hand, any dependence of the HMGs on the input PDFs and 
fragmentation functions is in principle artificial, since in the actual global analyses of SIDIS 
data the distributions would be determined uniquely and self-consistently by implementing 
the hadron mass corrected expressions in Eqs. (j38ll directly. 


23 


















IV. CONCLUSION 


With the imminent completion of the 12 GeV upgrade of Jefferson Lab, and with ongoing 
programs at existing facilities, a new generation of SIDIS experiments will vastly improve 
our understanding of the spin and flavor decomposition of parton distributions in the nu¬ 
cleon, and explore the relatively new domain of transverse momentum dependent parton 
distributions. A full utilization of the new data will require that effects associated with 
kinematical constraints at hnite energy are properly taken into account. Following earlier 
work which studied the dependence of unpolarized SIDIS cross sections on the masses of 
hadrons in the initial and hnal states 4ll-l43|. in this work we have presented a compre¬ 


hensive analysis of hadron mass corrections to both spin-averaged and spin-dependent cross 
sections and asymmetries at hnite values of Q^. 

Using the framework of collinear factorization, we have derived formulas for SIDIS cross 
sections in the presence of HMCs, which at leading order in as result in a rescaling of the 
PDFs in terms of the modihed Nachtmann variable and of the fragmentation functions in 
terms of the £nite-Q^ fragmentation variable (h- Our results respect all kinematical limits 
at hnite Q^, and reproduce the standard expressions in the Bjorken limit. An interesting 
feature of the modihed formulas is that, in the presence of HMCs, the parametric dependence 
on the scattering and fragmentation variables in the parton distribution and fragmentation 
functions becomes kinematically intertwined. While this complicates the analysis of SIDIS 
data in certain regions of kinematics, the ehects are calculable within our framework. 

We have quantihed the hadron mass ehects numerically as a function of the kinematic 
variables xb, Zh and in order to determine the regions where the corrections are most 
relevant. Generally, the HMCs are strongest at large xb and low (as for target mass 
corrections in inclusive DIS), and for large as well as very low values of the fragmentation 
variable Zh- The ehects are also more dramatic for heavier hadrons such as kaons than 
for pions at the same kinematics. Extending the previous analysis of mass corrections in 
unpolarized SIDIS j^|, we have also examined the dependence of the HMCs on the choice 
of scattered parton virtuality, In some cases the uncertainty due to this choice is quite 
signihcant, particularly at small values of and low Zh for pions, and over a larger Zh range 
for kaons, where the correction is observed to change sign. 

To illustrate the importance of HMCs in practical applications, we have computed the 
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corrections to SIDIS cross sections and polarization asymmetries that would need to be 


applied at kinematics relevant to upcoming experiments at Jefferson Lab |6lH63|. Here the 
xb and are necessarily correlated, so that usually the data bins at small xb correspond 
to lower values, while at large xb the is typically higher. For unpolarized pion 
production, the HMCs are strongest at large Zh, for all kinematics. At intermediate Zh 
values the corrections at low xb are relatively small, < 10%, but increase to ~ 40% —50% at 
higher xb {xb = 0.725), even at moderately large ~ 9 GeV^). Qualitatively similar 

behavior is observed for the semi-inclusive polarization asymmetry A^, which receives larger 
HMCs at higher xb values, although there is stronger sensitivity to the specific behavior of 
the input PDFs. Overall, our analysis suggests that mass corrections may be an important 
ingredient in future analysis of SIDIS data from facilities such as Jefferson Lab, especially 
at high values of xb, and particularly for hadrons heavier than the pion. 

An immediate application of the results derived here will be in upcoming global spin 
PDF analyses, such as by the JAM Collaboration which aims to £t an expanded set 
of high-energy scattering data, including SIDIS, down to = 1 GeV^. Future theoretical 
development of this work should include extending the calculation to next-to-leading order 
in as, which will necessitate consideration of hadron production at nonzero transverse mo¬ 
menta, ph±. The work can also be extended to other types of distributions measured in 
various single-spin asymmetries in SIDIS reactions, which provide information not only on 
the Xb and Zh distributions but also on the transverse momentum of the partons. 
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